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We study the entanglement in the ground state of a chain of free spinless fermions with a single 
side-coupled impurity. We find a logarithmic scaling for the entanglement entropy of a segment 
neighboring the impurity. The prefactor of the logarithm varies continuously and contains an impu- 
rity contribution described by a one-parameter function, while the contribution of the unmodified 
boundary enters additively. The coefficient is found explicitly by pointing out similarities with other 
models involving interface defects. The proposed formula gives excellent agreement with our nu- 
merical data. If the segment has an open boundary, one finds a rapidly oscillating subleading term 
in the entropy that persists in the limit of large block sizes. The particle number fluctuation inside 
the subsystem is also reported. It is analogous with the expression for the entropy scaling, however 
with a simpler functional form for the coefficient. 



I. INTRODUCTION 

The entanglement properties of many-body systems 
have attracted considerable attention and have become 
the topic of a large number of studies in the last 
decadeii^. In particular, the study of a suitable mea- 
sure of entanglement such as the entanglement entropy 
was triggered by the need to understand ground state 
properties that lead to the emergence of area laws^. This 
keyword refers to a rather generic property of ground 
states of bipartitioned systems in which the entropy S of 
a subsystem scales as the number of contact points with 
the environment. The most notable and well understood 
counterexamples are represented by one-dimensional crit- 
ical quantum systems with an underlying conformal sym- 
metry where the area law is violated in the form of uni- 
versal terms, scaling as the logarithm of the subsystem 
size^. In the translationally invariant case this can be 
written as 



S ■ 



■ In L + fco 



(1) 



with c being the central charge of the conformal field 
theory (CFT) and fco a non-universal constant. 

The presence of impurities has interesting effects on 
the entanglement entropy^. An example is given by crit- 
ical lattice models where single defective links separate 
the two subsystems. This can lead to a modified prefac- 
tor for the logarithm varying continuously with the de- 
fect strength in models with free fermions^"—, while for 
interacting electrons the defect either renormalizes to a 
cut or to the homogeneous value in the L —^ oo scaling 
hmitfl*''^^. In the more general context of a conformal in- 
terface separating two CFTs the effective central charge 
has been calculated recently and was shown to depend 
on a single parameter—. 

Another broad and intensively studied class of impu- 
rity problems is related to the Kondo effect^. From the 



viewpoint of block entropies a spin-chain version of the 
Kondo model has been studie d "'^^i^^ using density matrix 
renormalization group (DMRG) vaethods^^. Here a 
different geometry was considered with an impurity spin 
coupled to one end of a finite chain. The induced change 
in the entropy varies between zero and In 2 and its quali- 
tative behavior for various system sizes and coupling val- 
ues can be described in terms of an impurity valence bond 
picture. 

A counterpart of the Kondo effect can be found in a 
single impurity model introduced by Andersoi*^^. A sim- 
ilar, exactly solvable model was studied independently 
by Fanoi^ where the interaction of a discrete state with 
a continuum of propagation modes leads to scattering 
resonances. The exact form of these resonances is con- 
trolled by the couplings between the modes. Here we con- 
sider a simple geometry where an impurity is side-coupled 
to a linear chain of electron conduction sites. This set- 
ting can also be realized experimentally with gated semi- 
conductor heterostructure quantum dots. In accordance 
with theoretical predictionsSS, Fano resonances were de- 
tected as dips in the conductance measurements at low 
temperatures-^-. Despite the large amount of theoret- 
ical and experimental work on Fano resonances^, the 
question how they affect entanglement properties is still 
unanswered. 

In the present paper we address this question by inves- 
tigating the simplest model capturing the main features 
of the Fano resonance. The Fano- Anderson modeP^ is 
described by a free fermion Hamiltonian, thus standard 
techniques are available for the study of its entanglement 
properties^*. We find that the entropy scaling of a block 
neighboring the impurity is of the form ^ with a prefac- 
tor Coff which decreases monotonously as the parameters 
arc tuned towards the Fano resonance and depends only 
on a well-defined scattering amplitude. The numerical 
analysis leads to the same functional form of CeS that 
has recently been derived for simpler fermionic models 
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FIG. 1: (color online) Model geometry for (a) the infinite case 
and (b) the semi-infinite case. 



with interface defects^, and also seems to be closely re- 
lated to the one found for conformal interfaces^. If the 
subsystem contains an open boundary we find a rapidly 
oscillating subleading term in the entropy that persists 
even in the large L limit. 

In the following we first introduce the model and the 
geometries considered, along with the methods used to 
calculate the entropy. In Section IIIII we treat an impu- 
rity attached to an infinite hopping model, determine the 
correlation matrix and present results on the scaling of 
the entropy as well as the particle number fluctuations. 
The effects of an open boundary will be detailed in Sec. 
IIVI followed by our concluding remarks in Sec. [V] The 
derivation of an asymptotic form of the correlation func- 
tions and some of the necessary formulae for the semi- 
infinite case are presented in two Appendices. 



• in the semi-infinite geometry one has N sites with 
n = 1 , . . . , iV while uq is arbitrary and we have 
open boundary conditions 

The two geometries are depicted in Fig. [1] In the fol- 
lowing we set t = 1. After a Fourier transformation one 
has 



H 



y^ £aclca 



^q^q'^q 



+ y"tg{cld + d'<Cg) + eadU (3) 



wher 



wavenumbers are tq — 



eo — cosg. The couplings and the allowed 



infinite and tq = —gJj^-^simiQq with q = -jv+T" ^"^^ 
the semi-infinite case, respectively, where n runs over the 
corresponding site-indices. The Fermi-level €qp is set to 
zero by applying a potential eo = cos qp ■ 

The Hamiltonian (|3|) is known as the Fano- Anderson 
model and is diagonalized by introducing new fermionic 
operators /„ by the transformation^'^ 



c? = X! ^'i^t ' '^9 = X! ^^11'f' 



'Jq' 



(4) 



where 



'7-(e«) 



??±(e?) = e<? -Cd - S(eg±i5) (5) 



tqVq' 



Hqq' — Oqqi 

Cq - cq. 

with the self-energy term defined as 



iS 



(6) 



II. MODEL AND METHODS 

We consider a model of non-interacting spinless elec- 
trons described by a ID tight-binding chain that is cou- 
pled to an impurity. The impurity is represented by an 
additional site from which tunneling events can only oc- 
cur to a specific site no of the chain. The Hamiltonian 



n n 

- g {cl d + d^Cno) + edd^d 



(2) 



where c„ and d are the fermion annihilation operators 
along the chain and at the impurity site, respectively, 
while t and g are the corresponding tunneling strengths. 
The potential eo sets the filling of the chain and e^ is the 
site-energy of the impurity. 

We will consider two different geometries: 

• in the infinite geometry the A'^ -|- 1 sites are indexed 
as n = —A/2, . . . , N/2 while the impurity is located 



at the center of the chain (no 
periodic boundary conditions 



0) and we take 



E{eq±zS) = Y. 



*l' 



tq - tqi ± ib 



(7) 



The infinitesimal terms ±i(5 are needed to regularize the 
sum. Note, that the number of different q values equals 
the number of all sites, including the impurity. In the 
limit A — >■ cxD one has a continuum of the unmodified 
single-particle spectrum e^ of the chain while additional 
bound states can emerge as real solutions e of the equa- 
tion 



e-ed-S(e) = 0. 



(8) 



Our aim is to calculate the entanglement entropy S — 
— Tr(p In p) of a block of L sites neighboring the impurity, 
where p denotes the corresponding reduced density ma- 
trix. Since one is dealing with free fermions, this can be 
obtained through the eigenvalues Qi of the correlation ma- 
trix Cm, 



\ *^^n Cri 



restricted to the block 1 < m,n < L 



S 



^01nO-E(l-C01n(l-0) 



(9) 



III. INFINITE GEOMETRY 

We first consider the infinite geometry and calculate 
the matrix Cmn analytically; this result is in turn used 
to numerically obtain the entropy. The particle num- 
ber fluctuations are also considered where explicit calcu- 
lations are performed using an asymptotic form of the 
correlations. 



Correlation functions 



The correlations {c\^^Cn) are calculated by first going 
over to Fourier modes Cq and subsequently applying the 
transformation ([4]). The expectation values to be evalu- 
ated are then trivial and read {flfq') = Sgg'xiQ) where 
in the ground state xil) = 1 for the modes e, < and 
zero otherwise. Moreover, in the infinite geometry it was 
shown that Eq. ([8]) always gives two real solutions with 
e+ > 1 and e_ < — 1 corresponding to bound states above 
and below the band^^. The energies e± are obtained nu- 
merically by solving the quartic equation ^ with 



S(e±) 



± 



(10) 



These modes f± have to be dealt with separately in the 
transformation (|3]) with the factors 



i^± 



fJ'q± 



"e± - Erf 



(11) 



where N± is set by the normalization condition |i^±P -I- 
J2q l/^qiP — 1- Since (fLf-) = 1, one has a contribution 
C^„ = C^{m + n) to the correlations coming from the 
lower bound state. Setting tq = —g/\/N + 1 and taking 
N —> oo the sums are replaced with integrals over the full 
band [— tt, tt] and can be rewritten as contour integrals 
over the complex plane. Evaluating the residues one has 



C\l) 



2ei-e^ed-l 



-;/«- 



(12) 



with I = m + n for m,n > while the correlation length 
is defined by ^Z^ = arcosh(— e_). 

There are two contributions from the conduction band. 
First we have the translationally invariant terms 



C^„ = C°{m - n) 



sn\qF{m — n) 
■Kim — n) 



(13) 



that give the correlations without the impurity arising 
from the Sqq' term in ([B]). The remaining term gives rise 
to the correlation contributions C^„ = C^ {m + n) ; these 
can be evaluated by applying the same methods used 
for determining C^(l). The calculation yields the simple 
form 



C\l)^ 



'' —smSiq)smilq + Siq)) (14) 



where the scattering phases for g > are defined as 



tan 5(g) 



9 



Cd iiinq{€d-<iq)' 



(15) 



Here Yiq = ImE(eg -I- i5) is the imaginary part of the 
retarded self-energy with the real part being zero. 

Collecting all the contributions, the correlation matrix 
is given by 

The integrals in (|14p must be evaluated numerically and 
are shown on Fig. [Altogether with the contribution of the 
bound state for some fixed values of the parameters g and 
ed- The values of C^{1) are found to oscillate around the 
exponentially decaying curve of C^{1). Since they appear 
with opposite signs in (|16p . the average value will cancel, 
corresponding to the screening of the impurity induced 
localized state by the conduction electrons. 
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FIG. 2: (color online) Contributions to the correlations from 
the bound state C''{1) and from the conduction band C^{1) for 
the parameter values g = 0.2 and ea — —0.1. The inset shows 
the difference between these two contributions, compared to 
the asymptotic form H17|) represented by the dashed lines. 

The overall contribution is shown in the inset of the 
Figure. A careful analysis shows that the asymptotic 
behavior has the form of Friedel oscillations 

(l) = —-sinSFCoslqpl + Sf) (17) 
nl 

where only the scattering phase Sp = S^qp) at the Fermi- 
level enters. C^{1) is depicted by the dashed lines on the 
inset of Fig. [5] and shows a good agreement with the nu- 
merical data. The derivation of Eq. (J17p is summarized 
in Appendix A. It relies on an appropriate transforma- 
tion of the integrals in (fH)) and is valid for I ^ Iq where 
the length scale Iq is given in Eq. (jA7|) . 



B. Entanglement entropy 

The elements of the correlation matrix (TTHl) will be 
used to obtain the entanglement entropy of a block nu- 



merically as described in Section|TTl In some simple cases 
one can already give an answer by looking at the limiting 
form of the correlations. Taking Cd — > ±00, the impu- 
rity site will either be completely empty or completely 
occupied and therefore it becomes decoupled from the 
rest of the chain. One has 5{q) = while the contribu- 
tions from the bound state also vanish C^{1) — 0, hence 
Cmn = C^im — n) and the entropy will just be that of 
a homogeneous hopping model given by (P) with c = 1. 
Obviously, the limit of a vanishing coupling yields the 
same behavior. 

On the other hand, one could take the limit .g ^ 1 for 
very strong coupling. This corresponds to \S{q)\ -> 7r/2, 
that is, one has resonant scattering at every wavelength. 
The ground state will correspond to a singlet formed by 
the impurity and site zero that is otherwise decoupled 
from the rest of the system. The bound state correla- 
tions contribute only on site zero C^{1) — Sifl/2 while for 
sites ?7i, ri > one has Cmn = C^{m ~ n) ~ C^{ra + n), 
which is the form for a semi-infinite chain. The entropy 
of this geometry is known to scale logarithmically with a 
coefficient 1/6, which is half of the value that appears in 
the homogeneous case^i. 



S(L) 




FIG. 3: (color online) Crossover in the entropy scaling in the 
resonant case e^ = at half filling and for various coupling 
strengths g. The upper and lower-most dashed lines have 
slopes 1/3 and 1/6, respectively, acting as guides to the eye. 

The previous analysis shows that resonant scattering 
tends to decrease the coefficient of the entropy scaling. 
However, in general a Fano resonance is concentrated 
around a single wavenumber. This is expected to have 
the largest effect on the asymptotic behavior if the lo- 
cation of the resonance coincides with the Fermi-level 
I'^('?f)| — > 7r/2 yielding the condition e^ — > e^^ = 0. We 
take this value of the site energy in Figure [31 in which 
we plot the entropy for a range of values for the coupling 
g. We see that the entropy crosses over between the two 
limiting behaviors, giving a slope of 1/3 for small L and 
1/6 for large L when plotted against InL. Therefore, 
taking the block size large enough the impurity behaves 
like a cut in the chain. 



This effect is reminiscent of the Kondo screening mech- 
anism, where conduction electrons form a cloud around 
a magnetic impurity. The Kondo effect has an associated 
length scale^S- that being the distance beyond which the 
impurity appears to be screened off. In the present case a 
natural length scale also enters, which marks the asymp- 
totic regime of the correlations and could be expected to 
appear in the entropy scaling. From Eq. (|A7p one has 
lo ^ 1/g'^ for e^ = and qp = 7r/2, which seems to be 
consistent with the location of the crossover moving to- 
ward larger L in Fig. [3] as the coupling g decreases. In 
order to reliably verify the scaling of the crossover length 
one would have to consider larger blocks. Unfortunately, 
the achievable sizes were limited to L « 400 by the in- 
creasing numerical difficulty in evaluating the oscillatory 
integrands for the matrix elements (1141) . 

For L ^ Iq the screening cloud is effectively in a sin- 
glet state with the impurity electron and cuts the system 
in two parts. The coefficient of the entropy therefore 
renormalizes to the value 1/6 corresponding to a semi- 
infinite chain with an open boundary. Note, that a sim- 
ilar renormalization behavior was found for interacting 
electrons in the Luttinger liquid regime bisected by a 
hopping defect^iii. 

Our further numerical analysis shows that, for arbi- 
trary parameter values and i 3> 1, the entropy can be 
written in the form 



S{L) 



CcS 



InL 



(18) 



where the effective central charge CeS varies continuously 
between the values 1/2 and 1. The same behavior was 
found for free fermionic models with hopping defects^i^ 
where the dependence of CcS on the defect strength was 
determined by data fits. Recently, Sakai and Satoh de- 
rived the same form as Eq. p^ for the case of two confor- 
mal field theories coupled through a conformal interface. 
In their case they found Ccff in a closed form with depen- 
dence only on a single scattering amplitudei^. These con- 
formal interfaces describe a discontinuity in the compact- 
ification radii of two bosonic CFTs, which correspond to 
Luttinger liquids with unequal interaction parameters on 
the left and right hand side^^. 

Although the result fronJ^ is not directly applicable, 
a closely related form was recently found to describe the 
case of free fermions with simple interface defects^. We 
therefore attempted to generalize these results to the 
present model. Motivated by Eq. (ITT)) , we argue that the 
asymptotic behavior of S{L) and thus Coff should depend 
only on the scattering phase at the Fermi-level and a suit- 
able scattering amplitude can be defined as s = cosSp- 
Hence, we propose 



Coff 



+ — / u ('Vl + (s/sinhu)2 - l\ du (19) 



where 1/2 comes from the unmodified boundary while the 
integral is, up to a sign, the same as ir*iS, and describes 
the contribution of the impurity. 




FIG. 4: (color online) Effective central charge, obtained from 
fits of the numerical data, as a function of td and for various 
values of g and qp- The lines show the analytical form Eq. 
(|19[) evaluated for the corresponding values of the parameter 



matrix of a 2D strip with a defect line. In turn, an al- 
ternative representation of the integral in dT5|) via dilog- 
arithm functions was obtained. In our present model, 
however, a direct calculation of Cch through the single- 
particle eigenvalues seems to be a far more challenging 
problem. 




Figure 2] shows the resulting CcS obtained by fitting 
(fT8)) on different S{L) data sets together with (fT9|) evalu- 
ated as a function of e^ for various fixed values of the cou- 
pling g and filling qp. One has excellent agreement with 
the conjectured analytical form. The only visible devia- 
tions arise for e^ ~ which can be understood through 
the crossover phenomenon shown on Fig. [31 In this pa- 
rameter regime the asymptotic behavior sets in only for 
larger L and one has considerable finite-size corrections. 

The integral (IT51) depends on the model parameters 
only through the variable 



e\ sin^ qp 
£2 sin^ qp + g^ 



(20) 



Interestingly, the above formula exactly coincides 
with the transmission coefficient of the Fano-Anderson 
model2S with the Fermi-energy set equal to zero. There- 
fore, the effective central charge is directly linked to 
a simple physical quantity that, in experimental real- 
izations, is obtainable via conductance measurements 
at very low temperatures. Although the description 
of Fano resonances measured in various nanostructure 
experiments^^ typically require more detailed impurity 
models, the simple relation between Cog and the trans- 
mission coefficient might still survive in some parameter 
regime. This could open up the possibility of an indirect 
measurement of the leading term in the entropy. 

It should be pointed out that the impurity integral in 
Ccff has been obtained analytically for a transverse Ising 
chain in a DMRG geometry^. The calculation relies on a 
correspondence between transfer matrix spectra of classi- 
cal 2D Ising models and the single-particle eigenvalues of 
T-L in the reduced density matrix p ^ e~^ of the quantum 
chain. This analogy is well-known in the homogeneous 
case2i^ and can be extended by considering the transfer 



FIG. 5: (color online) Subleading term k of the entropy as a 
function of td for different fillings and g = 0.5. The homoge- 
neous value ko has been subtracted for better comparison. 

We also extracted the subleading contribution k in p^ 
from our data as shown in Fig. [5]for different values of the 
filling. Its value ko without the impurity also depends on 
qp^^ and was subtracted in the figure. One has a peaked 
structure around e^ = 0, however, unlike Ccff the constant 
k is in general not symmetric with respect to ed- Instead, 
one has the relation 



S{g, Ed, qp) = S{g, -e<i, tt - qp) 



(21) 



which is a direct consequence of the symmetry property 
Eq. (|A8|) of the correlations proven in Appendix A. 



Particle number fluctuation 



In the last part of this Section we will investigate the 
fluctuations in the number of electrons N = X]n=i '-n'^" 
contained in the block next to the impurity. Although 
being a simpler physical quantity, its scaling properties 
were shown to be similar to that of the entanglement en- 
tropy in a variety of ID quantum systems^ and, in the 
case of free fermions, also in arbitrary dimensionsMi^^. 
Therefore it is an interesting question whether this con- 
nection persists for the present impurity problem. 

The particle number fluctuation is a simple quadratic 
function of the correlation matrix and without the impu- 
rity is readily evaluated as'^^ 

(A7V2)o = trC"(l-C°)« ^InL + fc^ (22) 

TT 

with /cq = (In 2 -I- 7 -I- l)/7r^. Note that here the leading- 
order logarithmic behavior was seen to emerge from 



the large distance decay properties of the correlations. 
Therefore, we will use the approximation in Eq. (1171) to 
write C « C" + C^. Carrying out the traces, one finds 
that the additional term tr (C^)^ has a logarithmic con- 
tribution while trC-'^(l — 20^^) evaluates to a constant. 
The fluctuations then read 






k' 



(23) 



with the prefactor given by Keff = (1 + cos^ i5f)/2. 

The result is tested by comparing with the fitted values 
of Kcff as shown on Figure IHl The agreement is very 
good with the only sizeable deviations appearing in the 
region \5f\ ~ 7r/2 where the Fano-resonance is close to 
the Fermi-level and the approximation (1171) breaks down 
for the numerically attainable block sizes. Note, however, 
that in general the small distance deviations from the 
asymptotic form lead only to a different numerical value 
of k' as compared with the one obtained by evaluating 
the subleading terms in the traces. Finally, one has to 
point out that, apart from the similar logarithmic scaling, 
the prefactor Koff of the fluctuations is given by a much 
simpler function of the scattering amplitude than that 
governing the entropy. 
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FIG. 6: (color online) Fitted coefficients KcH of the logarith- 
mic term in the particle number fluctuation for various values 
oi g and qp, plotted against the scaling variable 5f- The dot- 
ted line is the function (l-fcos^ 5f)/2. The dashed line shows 
Cafi as a function of 5f for better comparison. 



IV. SEMI-INFINITE GEOMETRY 

We now turn to the investigation of the semi-infinite 
model. Our motivation behind studying a system with 
an open boundary is two-fold. On one hand, we would 
like to cross-check our argument used for explaining Eq. 
(|19p and show that the contributions in Coff are additive. 
In other words, we want to test our expectation Cog — 
Coff — 1/2 for the semi- infinite effective central charge Coff 
of a subsystem with an open end. On the other hand. 



boundaries were shown to induce interesting subleading 
behavior of the entropy for pure chaina^, therefore it is 
interesting to extend these investigations to the present 
impurity problem. 



A. Correlation functions 

The main difference from the infinite case is that the 
couplings tq are now dependent on the wavenumber. In 
principle the calculations are very similar, leading to 
slightly more lengthy expressions that are therefore sum- 
marized in Appendix B. However, there is an additional 
feature which leads to the emergence of bound states in 
the continuum (BIG) which were studied before^i^S. We 
will show how these BIG naturally enter the problem at 
the correlation function level. 

Our main interest is to determine the entropy of a block 
with L = riQ — l sites located on the left hand side of the 
impurity. Apart from the usual bound state contribution, 
the correlations from the conduction band now include 



^vn.n. -^ 



IF 



Aq 



TZ{q) sin qm sin qn (24) 



with m, n < riQ and we have defined the function 



n{q) = 



(e« - ed)' 



(eq-e,-Aqr+Tl 



(25) 



where A^ and Tq are the real and imaginary parts of 
the self-energy (|B5I) . respectively, and are related to the 
parameter Eg — ImE(eg + id) of the infinite system as 



Eg sm2noq 



E„ 2 sin noq ■ 



(26) 



Note that the term C^„ does not now appear in the 
correlations. 

It is instructive to analyze the main features of the 
integral in Eq. ((24)) . Taking \ed\ — >■ oo one has Tl{q) = 1 
and the integral reproduces the correlations of a pure 
semi-infinite chain. For very strong coupling g ^ oo 
one has Eg — >■ oo and therefore both of the parameters 
in ([25)) take very large values apart from a discrete set 
defined by the condition Qn = mr/no, n = 1, . . . ng where 
Ag = Fg = 0. Therefore, the function TZ{q) will become 
sharply peaked around these g„ values and formally one 
can substitute 



/^7^(,)^^5: 

J TT Uq^ 

9n 



(27) 



where the factor l/rio is needed for proper normaliza- 
tion. Setting no = L + I, the qn correspond exactly to 
the allowed wavenumbers in a finite chain of length L. 
Therefore the block is effectively decoupled by the im- 
purity and Crnn reproduces the correlations in a finite 
segment. 

The above argument breaks down for special g„-values 
fulfilling eg„ = ed where TZ{q) becomes zero and the cor- 
responding peak is missing. However, this is exactly the 



S(L) 




FIG. 7: (color online) Entropy scaling for a segment with an 
open boundary for different values of ed at half filling and 
g = 0.5. The dashed lines have corresponding slopes Ccff = 
Cefi — 1/2 and are shown for comparison. 



while r 



7r/2 



for no even and F, 



/2 



2E^/2 for no odd, 



and thus oscillates around the average E.;r/2- This effect 
could result in subleading corrections to the entropy that 
survive the no — > oo limit. 
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condition for the existence of the BIG. It has been shown 
that these states have non-vanishing amplitude only at 
the impurity site and inside the segment where they re- 
produce, up to normalization, the eigenstates of a chain 
of length L = riQ — 1^ . For 5 — > oo even the normaliza- 
tion becomes exact and the contribution of the missing 
peak is therefore reincluded this way. 

While for these special values of e^ the BIG are exactly 
located at g„, the resonances of Tl{q) gradually shift to- 
wards the wavenumbers of the decoupled block as the 
coupling to the impurity becomes larger. Hence, their 
contribution to the entropy is expected to diminish. 



B. Entanglement entropy 

The numerical study of the entropy is now carried out 
using the exact form of the correlations Cmn + C^„. The 
scaling of the entropy at half filling is shown on Fig. [7] 
for various e^ values. In agreement with our expectations, 
the slope of the curves coincides well with CeS = CeS — 1/2 
as illustrated by the dashed lines. Additionally, one ob- 
serves large amplitude oscillations in the data which seem 
to persist for large values of L. Although such oscilla- 
tions were already pointed out in case of quantum chains 
with a boundary, they were found to decay according to 
a power law and therefore vanish in the limit L — >■ 00^. 

A qualitative argument for this alternation can be 
given by considering the main parameters (|26p enter- 
ing the integrals. They contain oscillatory functions of 
q that, for large no, are expected to average out upon 



integration, yielding Aq = and F^ = Eq. These are 
just the values for the infinite case, explaining the aver- 
age behavior of the entropy. However, as seen before the 
Fermi level plays an important role in the asymptotics 
and at q = qp the parameters assume values which dif- 
fer from the average. For half-filling one has A7r/2 = 



FIG. 8: (color online) Alternating term AS in the entangle- 
ment entropy at half-filling, as obtained from data fits. The 
symbols correspond to different values of e^ and g and show 
a nice collapse when plotted against the scaling variable 5f 
of the infinite geometry. 

To extract this term, we have fitted our data according 
to S{L -t- 1) - S{L) = AS + a/L. This form, in general, 
gives a good description of the alternating part. The re- 
sults on AS' are plotted on Fig. [8] against the variable 
Sp- This choice is justified through our previous argu- 
ment, where the only non- vanishing parameter is given 
by 2E7r/2 and thus the scattering phase of the infinite 
case is still expected to be a relevant scaling variable. 
The data shows indeed a nice collapse to an interpolated 
scaling function which is shown by the dotted line. It 
has a zero a,t Sp = corresponding to the pure semi- 
infinite chain where the alternating term is known to 
vanish asymptotically^!. Note, that the symmetry un- 
der the exchange 5p — > —Sp is valid also for S{L) itself. 

The behavior of the entropy becomes more complicated 
when we choose another value of the filling. This is il- 
lustrated on Fig. [S] where we have compared qp = 7r/3, 
corresponding to one-third filling, with qp — 1 which 
gives an incommensurate value for the filling. In the for- 
mer case the entropy curve splits into three parts, shifted 
from each other in a similar way as was observed for half- 
filling on Fig. [T] This corresponds to the three possible 
values A^/3 and Ftt/s can take. However, in the incom- 
mensurate case the number of distinct values is infinite, 
resulting in a highly irregular entropy behavior. Since 
7r/3 ~ 1.05, the effective central charges for the two cases 
are almost equal which explains why the averaged slope 
of the curves is very similar. Nevertheless, for the small 
block sizes numerically achievable, the amplitude of the 
oscillations has the same magnitude as the average en- 
tropy. Therefore, a quantitative understanding of the 
subleading term would be clearly desirable and requires 
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further investigation. 
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FIG. 9: (color online) Oscillating behavior of the entropy for 
a commensurate qp — tv/S and an incommensurate qp ~ 1 
value of the filling. The parameters are g = 0.5 and ej — —0.8. 



The investigation of the semi-infinite geometry has, on 
one hand, supported our general behef that the leading 
contributions from both of the interfaces bordering the 
segment simply add up in the entropy. On the other 
hand, we found a subleading term which shows an inter- 
esting oscillatory behavior. The amplitude of the oscil- 
lations remains finite even asymptotically, an effect that 
has not been observed for pure chains. Although some 
arguments were given to describe its main features, a 
proper understanding could only be achieved through, 
analogous to the infinite case, a careful investigation of 
the asymptotic correlations. 

Finally, it would be worth extending this study to the 
case where the impurity is located further apart from the 
boundary of the subsystem. If this distance grows large, 
the effective central charge must return to the value of a 
pure system and an interesting crossover behavior might 
emerge. 
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V. CONCLUSIONS 

We have studied the entanglement in the ground state 
of a single impurity problem described by the Fano- 
Anderson model, focusing on the analysis of the effec- 
tive central charge that appears in the entropy scaling. 
We provided strong numerical evidence for the validity 
of a formula which gives the functional form of CeS for 
arbitrary parameter values of the model. In particular, 
this formula depends only on a single parameter, given 
by Eq. (j20p . that has a simple physical interpretation, 
namely it is the transmission coefficient of the impurity. 
This, in turn, establishes a connection between the entan- 
glement entropy and the low temperature conductance of 
the quantum chain that can be measured in experiments 
with quantum-dot nanostructures^^. Therefore, it would 
be interesting to generalize the study of the block entropy 
for Anderson-type impurities, including on-site Coulomb 
interaction between electrons of different spin, that give a 
more realistic model of the quantum-dots in experiments 
and yet are still tractable with DMRG techniques. 

On the other hand, an exact analytical treatment of 
the present problem is still lacking. One could follow the 
lines of Ref,^ and try to relate the reduced density ma- 
trix to the transfer matrix of an appropriate 2D classical 
system. Another possible approach would be the direct 
calculation of the entropy from the correlation matrix, 
which, in the pure case, requires knowledge of the eigen- 
values of a Toeplitz matrix^^ . With an impurity one has a 
sum of Toeplitz and Hankel matrices, which poses a much 
more difficult problem. However, the analysis of the cor- 
relations has led to a simple asymptotic form, enabling 
us to calculate the particle number fluctuation that in- 
volves, instead of diagonalization, only the traces of the 
correlation matrix. 
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Appendix A: Asymptotic form of the correlations 

The expression for the correlation functions involves 
the integrals in Eq. (J14p containing an oscillatory inte- 
grand. In the limit I 3> 1 one could approximate it by 
considering the scattering phase to be a constant. How- 
ever, this method only works for slowly varying S{q) func- 
tions which is not fulfilled by (J15p for g w and eq ^ ed- 
it will be shown, that with a suitable transformation of 
the integral one can extract the resonant contributions 
which exactly reproduce the term C''{1) associated with 
the bound state. To show this, we first rewrite ((T4| in 
the form 



C\l)= / ^ -.-"-^-^-.v^.-^i----'^- (Al) 
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(A2) 



Note, that the denominator is just the product 
r\j^{eq)r\-{eq) which is strictly positive along the integra- 
tion path in (jAll) but in general has four roots on the 
complex plane2&. Introducing the variable z = e*'' the in- 
tegration is now taken along an arc of the complex unit 
circle which can be written as the sum of the two con- 
tours shown in Fig. 1101 After changing variables and 




FIG. 10: (color online) Integration contours and poles of 
the integrand appearing in (|A3|) . The dashed line shows the 
boundary of the complex unit disk. 



which marks the regime I ^ Iq where the approxima- 
tion is valid. Note, that ^o can only grow large if the 
conditions e^ -^ g^ -^ 1 are satisfied. 

The contour integral can also be used to prove an ex- 
act symmetry property of the correlations. We define 
the particle-hole transformed functions C''{1) and C^{1) 
by performing the changes e^ -^ —Cd and qp ^ tt — qp- 
Now, one can show that the sum C^{1) + (— l)'C'^(Z) can 
be written as the integral (jA3[) over the complete unit cir- 
cle which is exactly evaluated as the sum of the residues 
at Z-|_ and z_. Using the symmetry property e+ = — £-, 
the latter pole contributes (— 1)'C'''(Z) which, after re- 
ordering, yields the relation 



C\l) - C\l) ^ (-1)' [C\l) - C\l)] 



(A8) 



using the symmetry under 2 — > z ^ it takes the following 
form 



c'{i) 



C1+C2 



dz ig'^z I 

2^i e(z)^ 



(A3) 



with the fourth order polynomial 

e(z) ^z'^ + 2z{7? - l)(ed - eo) + V^^ - 1 (A4) 

It is easy to show, that the polynomial has two real 
solutions z± = ±e~^/^T which lie inside the unit disk 
and 2+ (z_) is related to the lower (upper) bound state 
solutions as e^ = —{z± -t- z±^)/2 — ^cosh^^^. The 
remaining solutions form a complex-conjugate pair Z3 = 
z| and lie outside the unit disk. The closed contour C\ 



therefore only encircles the pole at z^ 
evaluated as 



and the residue is 



T, 4(7^Z , 

Res -^7} 

2=2:+ Kd\z) 



C\l) 



(A5) 



where C^{1) is defined under Eq. (jl2[) . The remaining in- 
tegral is taken along the contour C2 which is parametrized 
and using the analytic continuation of 



as z 



giig^g- 



the scattering phase it reads 



Re 



dq' 



smSiqp + ig')e''(«^+^«')'e''^(«^+'"'') (A6) 



Note, that S{qF + iq') varies smoothly and vanishes for 
g' — > cxD. For I 3> 1 the integrand is rapidly decaying, 
thus one can approximate the phase around q' = as 



S{qF + iq') ~ 5f -\- iq' 



J AS(q) 
dq 



. If the derivative is small. 



that is we are far away from the Fano resonance, one can 
set S{qp + iq') = 6p and the integral (jA6l) can be trivially 
carried out to yield the result in Eq. (|T7)) . The derivative 
term is then used to define a length scale 



d6iq) 



dq 



QF 



g'^jsir? qp - edcosqp) 
el sin^ qp + g'^ 
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Appendix B: Semi-infinite formulae 

We present here some of the formulae which is nec- 
essary for the calculations in the semi-infinite geometry. 
First, the self-energy function outside the band reads 



E(e±) = ± 



1 



|e±l 



-1 



2no 



and the bound state energies e+ > 1 and e_ < 
obtained numerically as the roots of 

e± - Ed- S(e±) = 0. 



(Bl) 



-1 are 



(B2) 



Note, that e+ (e_) exists only for parameter values ful- 
fiUing Erf > 1 - 2^^7X0 {cd < -1-1- 2g'^no). On the other 
hand, using the definition of the inverse correlation length 
e± ~ ±cosh^^ the second term in the parentheses in 
(jBip can be written as e^^"°/'»± . Therefore, if the impu- 
rity lies deep in the bulk of the chain with no — > 00 one 
has S(e±) = S(e±) and consequently the same solutions 
for e± as in the infinite case. 

The contribution C^^j from the lower bound state can 
be evaluated as 



g^N_ j4sinh2|^e-(™+")/«- 
e?.-l J4sinh^sinh^e-2"o/«- 

with the normalization factor 



iV_ = 



if TO, n > ftp 
if m, n < no 



2ei_ - e-ed - 1 - 2no52e-2»o/?- 



(B3) 



(B4) 



The limit no — > c» gives N^ = iV_ and the right hand 
side of (jB3|) becomes e"'-™ +" •'/'>" with the shifted indices 
m! ~ |m — no| and n' — |n — no|, thus recovering again 
the result P^ for the infinite case. 

Finally, one needs the expression for the retarded self- 
energy inside the band 



S(eo 
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which is an oscihatory function of the position no- This 
eventually leads to the emergence of bound states in 
the continuum, as described in the text. The scattering 
phases are defined as 



taiid{q) 



r„ 



Cd 



(B6) 



with A, = Re S(eq + iS) and T, = Im E(e, + iS). 

The band contributions also have to be treated sepa- 
rately on either sides of the impurity. On the right hand 



side, analogously to the infinite case, it can be expressed 
in the form C^(m — n) — C^{m + n) with 



C\l) 



-^ cos{ql + 25{q)) (B7) 



while the translationally invariant term C^{m — n) — 
C'^{Tn — n) is unmodified. The correlations on the left 
hand side cannot be given in that simple form and are 
analyzed in the text. 
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